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ABSTRACT 

In this paper, we propose a new adaptation of the D-iteration 
algorithm to numerically solve the differential equations. 
This problem can be reinterpreted in 2D or 3D (or higher di- 
mensions) as a limit of a diffusion process where the bound- 
ary or initial conditions are replaced by fluid catalysts. Pre- 
computing the diffusion process for an elementary catalyst 
case as a fundamental block of a class of differential equa- 
tions, we show that the computation efficiency can be greatly 
improved. The method can be applied on the class of prob- 
lems that can be addressed by the Gauss-Seidel iteration, 
based on the linear approximation of the differential equa- 
tions. 

Categories and Subject Descriptors 

G.1.3 [Mathematics of Computing]; Numerical Analy- 
sis — Numerical Linear Algebra 

General Terms 

Algorithms, Performance 

Keywords 

Numerical computation; Iteration; Linear operator; Dirich- 
let; Laplacian; Gauss-Seidel; Differential equation. 

1. INTRODUCTION 

The iterative methods to solve differential equations based 
on the linear approximation are very well studied approaches 
[S], [H, [IS, g], [IZ], [5], [IS]. The approach we propose 
here (D-iteration) is a new approach initially applied to nu- 
merically solve the eigenvector of the PageRank type equa- 
tion [11], [IQ], [9], [2], [8], [H]- 

The D-iteration, as diffusion based iteration, is an iter- 
ation method that can be understood as a column-vector 
based iteration as opposed to a row-vector based approach. 
Jacobi and Gauss-Seidel iterations are good examples of 
row-vector based iteration schemes. While our approach 
can be associated to the diffusion vision, the existing ones 
can be associated to the collection vision. 

In this paper, we are interested in the numerical solution 
for linear equation: 

A.X = B (f) 

where A and B are the matrix and vector associated to 
the linear approximation of differential equations with initial 
conditions or boundary conditions. 



In [6] , it has been shown how simple adaptations can make 
the diffusion approach an interesting candidate as an al- 
ternative iterative scheme to numerically solve differential 
equations. In this paper, we propose a new approach based 
on the pre- computation of the elementary diffusion limit. 
This limit can be then used for a given class of differential 
equations, for instance for 2D and 3D case, or for higher 
dimension. 

In Section [2] we introduce the 2D problem formulation. 
In Section [21 we define the notion of catalyst position and 
elementary solution. Section [4] describes the algorithm with 
the use of the elementary solution. Finally, Section [S] gives 
an illustration of the application and an evaluation of the 
run time gain. 

2. FROM THE HEAT EQUATION 

A typical linearized equation of the stationary heat equa- 
tion in 2D is of the form: 

r(n, rn) = ^ {T{n - 1, m) + T(n, m - 1) 

+T{n+l,m) + T{n,m + l)) 

which can be obtained by the discretization of the Laplacian 
operator in Cartesian coordinates: 

m^,y) = + (2) 

inside the surface (for instance, Q = [0,Lx] x [0,Ly]). 
Then additive terms appear for the initial or boundary con- 
ditions (Dirichlet) on the frontier dfl (for instance, for a; = 
or y = etc). 

More generally, we consider here linear equations of the 
form: 

T{n,m) = ^{T{n-l,m)+T{n,m-l) (3) 

+T{n+l,m) + T{n,m + l)) (4) 

when (n, m) E — Q \ and 

T(n,m) = g{n,m). (5) 

when (n, m) G dfl. For the general case, dfl is not necessar- 
ily only the boundary of Q (at least for discrete formulation, 
its asymptotic limit to the initial continuous problem is an- 
other story) and we may add any set of points included in 

n. 

We recall that the D-iteration requires updating two vec- 
tors (cf. [10]): the fluid vector F and the history vector H 



instead of a single vector for the Gauss-Seidel. The above 
equation can be solved or by iterating the equation (Q on 
or by applying the diffusion process associated to the 
D-iteration. 

3. ELEMENTARY SOLUTION 
3.1 Diffusion on ID 

For the purpose of illustration, let's consider the ID case 
with the following equations (from differential equation of 
order two): 



Then let's consider the limit of ([6]) when we have an el- 
ementary catalyst: by symmetry, we can just explore the 
space IN: 

• at the first iteration, the half of 1 is sent to the position 

1; 

• at the second iteration, half of 1/2 (1/4) is sent to 
positions 2 and 0: is a catalyst, therefore the 1/4 it 
receives disappears; 

• at the thirst iteration, 1/8 is sent to 1 and 3, and so 
on. 



T(n) = -(T(n-l)+T(n+l)) 



when n G Q," and 



T{n)=g{n). 



(6) 



(7) 



0. 



when n € dfl. 

The solution can be found by iteration of equation (|6]) on 
fi" with boundary condition: 

T{n) = g{n) when n G dQ 
T{n) = when n ^ dO.. 

This is exactly the Gauss-Seidel iteration if we always use 
the latest updated values of T and the limit is the piecewise 
linear function joining the points {n,g{n)) at n £ di}. 

The diffusion based iteration would apply the equations: 



H{n)+ = F(n); F{n) = 0, (8) 
F{n + 1)+ = ii^(n), F{n - l)+ = iF(n). (9) 



3.2 Catalyst position 

We first introduce the notion of catalyst node: a node 
(position n) is a catalyst if it diffuses once its initial fluid 
value F{n), then behaves as a black hole, i.e., it absorbs all 
fluids it receives without retransmitting them to its neigh- 
bour positions (so that we have H{n) is constant to F{n) 
after its diffusion). 

Theorem 1. If we associate dQ. to catalyst positions imth 
the boundary conditions ^ defining their initial fluid val- 
ues, the D-iteration diffusion's limit is equal to the limit of 
the iteration of the equation \^ 

Proof. This theorem is a direct consequence of the equa- 
tion on H which does exactly ^ (cf. |12)V 

To solve efficiently such an iteration scheme (of course, 
the interest is for 2D/3D or more complex mix of differen- 
tial equations of higher order in 2D/3D), we introduce the 
notion of elementary solution for the fluid diffusion process: 
the elementary solution is the limit of the diffusion process 
when we put a fluid 1 at position (0, 0) and with boundary 
condition zero on the frontier. If we impose that (0, 0) is a 
catalyst, the solution can be obtained from the elementary 
solution by re-normalizing all values by 1/(1 — H{Q, 1)). 

We call elementary catalyst when we have dO. = {0, A'^} 
and g(0) = 1, g{N) = with ^ = [0, iV]. 

We can compare this idea to the use of the Green's func- 
tion to solve linear differential equations. The difference is 
that the function we define here is much easier to compute 
and not dependent on the boundary shape. 
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Figure 1: Elementary catalyst on ID: diffusion on 
IN. 

If we iterates the diffusion, at the limit we find a discrete 
function which can be used to solve all equations of type ^ 
(see Section [4]). What's interesting is that when we don't 
put the boundary A*', we can prove that the amount of fluid 
that disappeared in the black hole is associated to the limit 
of hypergeometric series deflned by: 

_ 2n + 1 

and converges to 1/2 (use of Gamma function, [2]). Since 
1/2 is the fluid sent to the direction IN, that means that all 
fluid are flnally absorbed by the black hole and that each 
value H{n) is convergent (non-decreasing and bounded). 
What's even more interesting is that at the limit we have 
H{n) = 1 for all n G IN (proof by induction from what's 
received at 0), which means that all positions will receive 
and send exactly 1 fluid: this is a consequence of the fact 
that the associated random walk is recurrent (only for ID 
and 2D). However the convergence speed to the limit of this 
process (by diffusion or by collection) is very slow (for large 
n). This is why it may be interesting to pre-compute those 
iterations once, but not up to its limit (which is the constant 
function), but with a boundary condition at N. 

In the equation we consider, we'll always have the fron- 
tier of n that are all catalysts, therefore all fluid reaching 
the border of Q. will disappear and this guarantees a faster 
convergence compared to the case without boundary. Since 
the limit is the limit of the diffusion iterations when all cat- 
alysts have injected exactly gin) and when [-/^(n)] = 0, we 
can use the catalyst limit we pre-compute on a flnite set as 
a common block for faster diffusion (see Section U). 

Figure [2] shows the limit function we obtained after 10^ x 
2000 iterations (10^ cycles over [1, ..2000]): note that this is 
far from linear function! As mentioned above, the limit is in 
this case the piecewise linear function that joins (—2000, 0) 
to (0, 1) and (0, 1) to (2000, 0). 
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Figure 2: Diffusion of the ID elementary catalyst on 

[-2000,20001. 



From the linear algebra point of view, this means that 
we have a matrix associated to the equation (|6} with initial 
condition g{0) — 1 and that its spectral radius is strictly less 
than 1 thanks to the boundary condition. 

We can also interpret this limit (after normalization) as 
the average sojourn time of a random walker when starting 
from position and to which we forbid to return to or 
reach the boundary position TV. 

Remark that, if was not a catalyst and sends back the 
fluid it receives, we end up with unbounded quantities (for 
H) when the boundary is not set, since the initial fluid 1 
never disappears (null-recurrent random walk, stating from 
0). 

Finally, note that in the general case, we will have equa- 
tion of the form (instead of Equation (O): 



T{n) = aT{n - I) + l3T{n + 1). 



(10) 



An example is given in the next section. We may have also 
equations where a and /3 may depends on n. An example 
of this is the discretization of the Laplacian in polar coordi- 
nates, invariant by rotation: 



dr^ r dr 



(11) 



which gives: 



3.3 Application on ID 

Consider a differential equation of 2nd order: 

y"(x) + ay'{x) + l3y{x)^ fix) (13) 

with boundary condition y{0) and y[Lx)- The naive corre- 
sponding iteration scheme (from discretization) is: 



Vn 



2 -f ae - /3e2 

((1 + Qe)y 71-1-1 + Un-l 



Vn 



(2-/362) 

((l-hae/2)i/„+i-f (l-ae/2)y 



(14) 
(15) 

(16) 

^Vn) (17) 



for increment e (y„ = y(ne) and /„ = f{ne)). To solve this 
iteration scheme, we can: 

• solve the elementary catalyst pre-diffusion; 



• then use the catalyst pre-diffusion to: 

— for each position ?i, diffuse e^/n; 

— solve the diffusion problem with boundary y{0) 
and y{Lx). 

For the diffusion of e^/„ we can use a pre-diffusion as for 
elementary catalyst but without the constraint of the initial 
position behaving as a black hole. Then the elementary 
catalyst diffusion limit can be obtained by normalizing all 
values by -ff(O). 

3.3.1 Simple example 

Assume we want to solve the equation: 

y"{^) = fix) 

with /(x) = {~0.99cos(x) + 0.2sin(x)) exp""^''^'' and with 
boundary conditions: y{0) = l,y(50) — (the solution is 
y{x) = cos(x)e~"/i''). 

The discretized equation is (from Equation H16|) with a = 
/3 = 0): 



Vn =^(t/n + l -I- yn-l) - ^ fn 



(18) 



Since the catalyst limit of (|SJ is the piecewise linear func- 
tion, we can solve the equation as follows: 

• for each position n, diffuse e'^g„/2, which means adding 
the linear function (multiplied by ~ 1/(1 — (Lx — 
1)/Lx), because {Lx — 1)/Lx is the quantity that comes 
back to the diffusion initialization point and here there 
is no black hole behaviour): 

for (int i=l; i < Lx-1; ++!){ 

transit = - step*step/2 * g(step*i) * (Lx-1) ; 
for (int j=-i; j < Lx-i; 
int x = i+ j ; 

Y[x] += transit * (Lx-1 - abs(j)) / (Lx-1); 

> 

> 

(where step is e); 

• solve the diffusion problem with boundary y(0) and 
yiLt): 

transit = 1.0 - Y[0] ; 

for (int i=0; i < Lx; ++i){ 

Y[i] += transit * (Lx-1 - i) / (Lx-1); 

} 

transit = 0.0 - Y[Lx-l]; 
for (int i=0; i < Lx; ++i)-[ 
Y[i] += transit * i / (Lx-1); 

y 

For any e we obtain directly the limit as a superposition 
of fluids diffusion from /„ and the boundary conditions and 
there is no need to do any iterations thanks to the explicit 
form of the catalyst limit. The results are shown on Figure 

m 




Figure 3: Example of ID. 

The explicit theoretical formulation of the algorithm in 
this example is: 

y{aL^) = {1 ^ a)y{G) + ay{L^) (19) 

+ 2i# (2m - i + \N/L^a - i| - aN) /{LJNi) 

(20) 

which in the limit (A'' — > oo) is for Lx = 1; 

y{x) = {1 ~ x)yiO) + xyil) (21) 

1 

+ -/ {2xt-t-x+\x-t\)y"(t)dt. (22) 
2 Jo 

This formulation can be directly solved if we look for an 
expression of y{x) as a function of x,y{0),y{l) and of the 
integral of the form J u{t, x)y" {t)dt. We'll see how such a 
formulation can be generalized. In this case, the diffusion 
approach is equivalent to the usual discretization of the in- 
tegral in the equation (|2ip . 

3.4 Diffusion on 2D 

In this section, we consider linear equations associated to: 

r(n,m) = i (r(?i- l,m) +r(n,m- 1) (23) 
+T{n + l,m) + T{n,m + 1)) (24) 
when (n, m) £ fi" and 

T{n,m) = g{n,m). (25) 

when (n, m) £ dQ. 

For 2D, we consider as for ID, the diffusion limit of the 
elementary catalyst at position (0,0). Figure |4] shows the 
limit function we obtained. 

The computation of this limit on a large space is compu- 
tation costly. Using the rotation invariant polar Laplacian 
equation, we can in fact find the explicit solution, which is 
of the form: 

^ 'og('') 

This function has a singularity at (because the limit to 
the continuous case must be a density or a measure). An 



Figure 4: Limit of tiie 2D elementary catalyst on 

[-100,100] X [-100,100]. 

empirical interesting candidate to approximate the diffusion 
limit of the elementary catalyst in 2D is: 

where a has been evaluate from the explicit diffusion it- 
erations. T{n) is then such that T(0) = 1 by definition, 
r(l) = a and T{Lr) = 0. In fact, we found that it was 
better to use the iteration of Equation (|12|) which was close 
to the above close formula, except the tails. However, those 
limits are in polar coordinates, which introduce a bias when 
applied on the Cartesian coordinates and which we can not 
eliminate (because diffusion on grid is not rotation invari- 
ant!). In this paper (Section [5|, we used first the iteration 
of Equation (|12p (which is ID, so very fast), then from this 
we defined the starting point of the iteration on 2D, using 
also the symmetry of 2D (computation on 1/8 of the plane), 
which accelerated the full naive 2D scheme iteration by fac- 
tor 10-50. 

Another alternative is to apply the ideas of Section |4] dur- 
ing the pre-computation: after iterations on a smaller space 
(for instance. A''' = N/2), we can save the results H'^,F'^, 
then we can replace the elementary diffusion by directly 
copying the results of N iterations as a block. This is inter- 
esting to gain an order of precision quickly, exploiting the 
fact that i^" has fluids concentrated at the border. However, 
after one block copy operations, we find again a configura- 
tion where the fiuids are spread more uniformly. Optimizing 
the pre-computation phase is an independent problem which 
we don't analyse further here. 

4. ALGORITHM 

We consider the 2D problem Ay — f on Q with boundary 
condition g on dQ for illustration. The method should be 
easily extended to a much general linear operator associated 
to other differential equations. 

We assume the elementary catalyst's limit is pre-computed 
on a finite set [—Lx,Lx] x [—Ly,Ly] with boundary condi- 
tion g{0,0) — 1 and g{x,y) = if = Lx or \y\ = Ly. 
For the sake of simplicity, we will consider of the form 
[0, Lx] X [0, Ly] (if not, we can choose Lx the maximal x- 
distance between two points of Q and similarly for Ly). For 



the practical computation, we iterate the D-iteration until 
the remaining fluid \F\ is below the targeted error. Then, 
we store in a file the last states on H and F (in the following 
denoted H° and F°). 

Then we apply the following process: 

• load the above results H° and F°; 

• define a new variable [Li:] [Ly] and FfLa;] [Ly] (initial- 
ized to 0); 

• set the initial fluid F equal to / in fi"; 

• diffuse F on Q, (including dfl): here, diffusion means 
adding J?" and i^" on H and F respectively at trans- 
lated position (by x and y); 

• diffuse fluid g{x,y) — H[x][y] on d^l (choose positions 
where \g{x, y)—H[x][y] \ is the largest or above a certain 
threshold; here, diffusion means adding and on 
H and F respectively at translated position (by x and 

y)-- 

Diffusion of "g(x,y)-H[x] [y] " : 
for (int x=0; x < Lx; ++x)-[ 
for (int y=0; y < Ly; ++y)-[ 

if ( bound[x][y] ){// boundary position 
transit = g[x] [y] - H[x] [y] ; 
if ( abs (transit) > Thresh_ ){ 
for (int i=0; i < n_x; i++){ 
for (int j=0; j < n_y; j++){ 
H[i][j] += 
transit*HO[abs(i-x)] [abs(j-y)] ; 

F[i][j] += 
transit*FO [abs (i-x)] [abs(j-y)] ; 
} 

} 



• the previous step is repeated until the threshold is be- 
low the targeted error; 

• if required, we may also diffuse fluid F which are above 
a given threshold (we may also decide not to use F at 
all), because as far as we keep F and H, all operations 
are invertible in the sense that we can inject the surplus 
or the deficit fluid to make the exact convergence in 
any order. 

The numerical solution to our problem is then given by 
H. If H is exactly equal to the boundary condition g on dfl, 
then // on f2 is the exact limit. 

As for the ID case, we can express this approach by 
the projection method where the elementary catalyst limit 
serves as a unique base. It can be also understood as an ap- 
plication of calculus of variations or a Lagrangian approach. 
Let's call <j) the limit of the elementary catalyst (for instance 
on a square surface that's bigger than i}). Then, we can 



rewrite the algorithm under the form: 
y{x) ~ y{n,m) 



y{n,m) = -^^fii 



(26) 
(27) 



+ iy{xb) ~ a{xt) ~ —^g{{)(fii{xt,)\ 4>^^ix), 

(28) 

where x = ((5n, (5m) e Q, (regular grid of 5), (^i{x') the 
value of at point x when the origin is set at i and ^ = 
((>/ (1 — 0(0, 1)), and a{xi,) is term expressing all diffusion re- 
ceived from other boundary condition related diffusion. Our 
approach can be understood as an iterative approach to find 
the coefficients a{xi,). 

Its limit (if existence) for 5 — >■ can be formulated as: 



(y{xb) — a{xb))(j){xb — x)dxb (29) 



- [ Ay{tmx-t)dt (30) 



/S.y{t)(f){xb — t)(f)[xb — x)dtdxb (31) 



an Jn 



where the second term comes from the diffusion of fluid in- 
side fl and the two other from the correction for the bound- 
ary conditions. This formula assume in particular that we 
have a limit of when A'' goes to infinity. We can inter- 
pret (j>{n, m) as the probability for 2D random walk to reach 
(n, m) before touching the boundary starting from (0,0). 
When A'^ goes to infinity, the random walk tends to the 2D 
Brownian motion and 4'{x) in the continuous space is the 
probability that from (0, 0) we reach [x, x -\- dx] x [y,y -\- dy] 
before the boundary is touched. 

In a particular case when / = 0, they is a very nice theory 
of probability which shows that y is given by an explicit 
integration formula ( [TSl [Til [3]): 

y{x) = h{x)=^4g{B{T))] 

where B is the Brownian motion, T the stopping time when 
the boundary is touched. If the boundary is a sphere, we 
have a more explicit formula of the form: 



y{x) = h{x) 



Sd-i \x - y\ 



■;9{y)^d{dy). 



From the diffusion point of view, we can understand why 
with the sphere we can have a simpler formula: the diffusion 
from one point of the sphere to all others points of the sphere 
follows exactly the same process, meaning that in our ap- 
proach the terms ^(a::;,) can be eliminated if (f> is associated 
to this diffusion model. 

Our approach can be understood as an explicit practical 
solution, not only in presence of /, but also for a general op- 
erators (so not only harmonic functions) associated to the 
differential equations, using a specific choice of <j). When the 
diffusion operator is not symmetrical in the four directions, 
the very nice theory of harmonic function does no more ap- 
ply. However, the idea of exploiting the pre-diffusion (0) can 
be also compared to the use of the Green's function G{x, s) 
(when it is known!) and express the solution as: 



y{x) 



I 



G{x, s)f{s)ds 



But while this is an exact solution, the computation of the 
Green's function may be even more complex than solving 
directly by an iterative scheme in a general case. 

Note that our algorithm has no guarantee of convergence 
(on a{xi,)). We hope address this point in a future paper, if 
such a consideration is not already proposed in the past. 

4.1 Error estimate 

The distance to the limit can be estimated from 
r= Yl \n^M\+ E \9ix,y} - H[x][y]\. 

The first component of r is the residual fluid resulting from 
the diffusion by catalysts and the second component is the 
surplus or the deficit fluid that are injected to O. If r = 0, 
H is the exact limit of the problem. 

5. EVALUATION 

5.1 Convergence comparison 

For the evaluation purpose, we considered the following 
(too simple!) scenario: 

• SI: a 2D diffusion problem (/ = 0) with Lx ~ Ly 
and boundary condition on the border: g{x, y) = 100. 
The solution of this problem is obviously a constant 
function equal to 100 on every point of SI. 

The results are shown on Table [T] we used 2 Linux lap- 
top: Intel(R) Core(TM)2 CPU, U7600, 1.20GHz, cache size 
2048 KB (Linuxl, g + + - 4.4) and Intel(R) Core(TM) 15 
CPU, M560, 2.67GHz, cache size 3072 KB (Linux2, g + + - 
4.6). The pre-computation of the elementary catalyst on 
[—Lx, Lx] X [—Ly, Ly] has been done for a given target error 
(target, on the remaining fluid); the runtime for this pre- 
computation is given by pre-comp. The results have been 
saved in a simple ASCII file, its loading time is given by Init. 
We observed that the limitation of the error of our approach 
was about 10~^ (which means for g{x) = 100 a relative pre- 
cision of 10~^), resulting probably from the double precision 
(about 10"^^)) we have on F° (relatively to H°). Through, 
this school case, we just want to illustrate the potential of 
our approach. 

5.2 Stationary heat diffusion in 2D 

Let's consider a simple variant of SI: we set 

• S2: a very simple diffusion problem with Lx = Ly — 
2000 and boundary condition on the border: g{0,y) = 
100 and g(x, 0) = g{x, Ly) = g{Lx,y) = 0. 

Results are on Figure [1 [H [3 [TOl [11] and [H for the 
D-iteration, we use the pre-computation (j> that is generated 
in the previous section. We can see that with the naive 
iterative method, the convergence to the limit may be really 
slow when a large grid is considered. The result obtained in 
30s with our approach has in this case a better convergence 
than with 16 hours with Gauss-Seidel (the gain is reaching 
a factor 2000). But of course, the gain was obtained thanks 
to the previous pre-computation (p which took about 1 day. 
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Table 1: Comparison of computation cost. Pre- 
comp: pre-computation time of the elementary cata- 
lyst, target: target error on the remaining fluid for 
pre-computation. Init: initialization time, error: 
distance to the limit, error 2: maximum increment 
of the last iteration for GS, r for DI. Lx = 100, 200. 
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6. CONCLUSION 

In this paper we addressed a first analysis of the poten- 
tial of the D-iteration when applied in the context of the 
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0.4 


0.14 


0.09 


error2 




le-'^ 


260 


270 


time 


9500 


123700 


53 


500 


gain 


1 


1 


xl80 


x250 


Linux 2 










Pre-comp 


X 


X 


6900 


93000 


target 


X 


X 


10"=* 


10"^ 


Init 


X 


X 


3 


10 


error 


1.0 


1.0 


1.0 


0.7 


error2 


le-5 


le-6 


1500 


2000 


time 


2800 


50000 


5 


30 


gain 


1 


1 


x560 


xl600 


error 


0.1 


0.1 


0.11 


0.09 


error2 


le-^ 


le-^ 


140 


270 


time 


4200 


55000 


30 


150 


gain 


1 


1 


xl40 


x370 



Table 3: Comparison of computation cost: Lx = 
1000,2000. Pre-computation is only done on Linux2. 



numerical solving of differential equations. We showed that 
using the regularity of the diffusion process, we can exploit 
the idea of the pre-diffusion. The diffusion approach gives 
a new way of understanding the differential and integration 
associated operator iteration at a fundamental level and of- 
fers a great potential for a very fast numerical computation. 
Further exploitation of this will be addressed in a future 
paper. 
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Figure 5: Gauss-Seidel. Run time: 4 min. 
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Figure 6: Gauss-Seidel. Run time: 1 hour. " 




Figure 8: Gauss-Seidel. Run time: 4 hours. 
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Figure 12: D-iteration: 30 s. 



Figure 9: Gauss-Seidel. Run time: 16 hours. 



